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Abstract. A family of discontinuous symplectic maps on the cylinder is considered. This 
family arises naturally in nonsmooth Hamiltonian dynamics and in switched Hamiltonian sys- 
tems. The transformation depends on two parameters and is a canonical model for the study 
of bounded and unbounded behavior in discontinuous mappings due to nonlinear resonances. 
This paper provides a general description of the map and presents one case of an unbounded 
orbit. 



1. Introduction. 

Theory of small perturbations of completely integrable Hamiltonian systems has a long history 
that goes back to 19th century effort to explain stability of planets. The major breakthrough 
occurred in the late 1960s, when Kolmogor ov- Arnold- Moser (KAM) theory was created. 

The KAM theory states that under some non-degeneracy conditions, stable motion persists in 
a completely integrable Hamiltonian system under sufficiently small and smooth perturbation. 

For the original application, 3-body problem, the smoothness was not an issue as the grav- 
ity force is analytic, outside of a small set of singularities. However, further applications of 
KAM theory to stability problems in physics and engineering, do require limited smoothness 
assumptions and also weaker forms of the so-called twist (nonlinearity) conditions. 

The degree of smoothness of the perturbation has a crucial role in the theory. In his famous 
ICM lecture Kolmogorov gave an outline of the theory where he required analyticity. Shortly, V.I. 
Arnold proved Kolmogorov's statement, also under the assumption of analyticity. Independently, 
combining Kolmogorov's method with Nash smoothing technique, Moser proved a KAM type 
theorem requiring 333 derivatives. Subsequently, the smoothness requirement was reduced to 
single digits (C^) and several counterexamples have been found for lower regularity maps, see 
e.g. 0. 

Moser proved his theorem for the case of area-preserving monotone twists maps of the annulus. 
In this article we also restrict our attention to the representative case of twist maps on the plane, 
which corresponds to the periodically forced Hamiltonian systems with one degree of freedom. 

The above KAM countrexamples, that were constructed for the general twist maps, do not 
provide a tool to decide stability in specific physics problems. Therefore, it is important to 
investigate special maps arising in physics and geometry. 

We note that even in the most extreme case of discontinuous maps, the stability problem 
is already nontrivial. In the next section, we review several such systems where boundedness 
problem for discontinuous maps naturally arises. Then we introduce a family of discontinuous 
twist maps, which captures the essential properties of those examples. The family contains a 
natural physical system which we call pinball transformation. The hallmark of the pinball map 
is the small twist, which on the one hand frequently occurs in applications, and on the other 
hand makes stability problem rather delicate. 



Date: March 25, 2013. 



2 



MAXIM ARNOLD, VADIM ZHARNITSKY 



2. Discontinuous twist maps and az-TRANSFORMATiON. 



Discontinuous maps arise naturally in Hamiltonian systems with impacts, such as Fermi-Ulam 
problem, billiards, and more recently in hybrid or switched systems. It is worth noticing that 
in each of these examples, under the additional smoothness assumptions, KAM theory applies 
assuring boundedness of energy. 

One should also keep in mind that while the general monotone twist maps are characterized 
by a function of two variables h{xi,X2), these particular examples correspond to symplectic 
maps characterized by function of one variable, e.g. for billiards h{xi,X2) = \\x2 — xi\\. Such a 
restriction makes it nontrivial to construct physically meaningful escaping trajectories. 

For the readers' convenience, now we briefly describe these systems. 

Example 1: Particle in square wave switching potential. Hybrid or switching systems is an active 
area of research in applied mathematics and engineering sciences, see e.g. [S])PP- A prototype 
example of a switching system, where boundedness problem is already non-trivial, is a classical 
particle in square wave periodic potential which changes the sign, periodically in time. 

More precisely, let the potential be V{x) = (— l)'^! and assume the potential is switched every 
second V{x.,t) = (—1)'*^ • (— 1)^^'. While such potential is not differentiable, there is a natural 
way to define the dynamics by using the energy relation: the kinetic energy changes by 2 if the 
particle passes f € Z integer points. It is common to ignore the singular subset of the extended 
phase space (i, i, x) where there is discontinuity in both time and space and the dynamics is not 
defined. Such subset has zero measure. Outside the singular set particle moves with constant 
speed V = ^JE ^1, gaining or losing energy by one at each switching. See the appendix for more 



Example 2: Fermi-Ulam accelerator. The Fermi-Ulam system consists of a classical particle 
bouncing between two periodically moving walls. The application of KAM theory shows that 
velocity (or energy) of the particle is uniformly bounded \v\ < C{v{0)), provided the periodically 
moving wall's position is sufficiently smooth p(t) G C^(0,T) [7J. 

Fermi-Ulam problem can be reduced to a particle traveling in a periodic non-smooth potential 



It turns out that lack of smoothness in x {e.g. due to the presence of the wall in Fermi-Ulam 
problem) does not destroy bounded behavior as one can exchange the role of time and coordinate 
and then obtain a smooth monotone twist map by integrating over x, see e.g. [8]. 

If there is lack of smoothness in both space and time in the periodic potential problem, then 
KAM theorems do not apply. 

In the worst case the map is discontinuous, but even then, finding unbounded solutions could 
be challenging. One case, however, is more tractable: if jumps in the velocity are so large that 
the solution makes full revolution over one period of forcing so it will be in tune for the next 
velocity increase. A typical example would be given by this map 



Such scenario takes place in Fermi-Ulam problem if p{t) has saw-tooth like shape. But, if 
the velocity increments are smaller, then the twist will eventually detune the solution out of 
resonance. 

Example 3. Outer Billiards. The question of boundedness becomes a lot more delicate and there 
are few examples of escaping trajectories for such systems in the dual billiards. Only recently, 
Schwartz and then Dolgopyat and Fayad constructed unbounded solutions in the presence of 



details. 



X + V'{x,t) = 0, where V' G L\T^). 




xi = {x + y) mod 1 
2/1 = y + sgn(xi - \). 
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piecewise smooth boundary. In the appendix, we give some heuristic description how the dis- 
continuous twist map arises in this example. 

az Map: A model of boundedness problem for discontinuous twist maps. 

In this paper, we introduce a two-parameter family of discontinuous monotone twist maps that 
seems to capture the essential difficulties of several switching-like (discontinuous) systems. The 
map is given by 

I xi = (x + ay^) mod 1 

\yi = y + sgn(xi - I) 

and will be referred to as az-map. Notice, that the map is invariant with respect to the scaling 
that varying the amplitude of the changes in the second variable or varying the length of the base 
circle in the first variable will led to the equivalent system with different values of parameter 
a. We also observe that az transformation preserves the unit-step lattice in action variable. In 
other words, the action variable is quantized for any fixed initial condition. 

For different values of parameters az map corresponds to some natural systems: 

• z = 1, Fermi-Ulam with saw-tooth p{t), discontinuous standard map. 

• z = 0, Kesten system (skew product of irrational rotation with jumps). It is defined 
in the next section. 

• z = 1/2, Particle in switching square wave periodic potential 

• z = -1, Pinball problem (studied in this paper) 

We explain in more details how az-transformation arises in each of these examples in the 
appendix. 

Zero twist example. Kesten system. The following system was introduced by Kesten in j6] 
independently of any connection with KAM theory- type of problems. He considered irrational 
rotation on the circle and asked what is the discrepancy between the orbit visiting different open 
subsets of the circle having equal measure. In particular one can consider different halves of the 
circle x £ (0, 1/2) and x G (1/2, 1). In our notation, his system corresponds to the map with 
z = 0. 

In this degenerate case, there is no twist in the system and the dynamics is a skew product. 
Thus, one can easily provide a set of values of parameter a (e.g. a = 1) for which there are 

unbounded orbits. On the other hand for a = - any trajectory of the system Q is bounded 

since any point has period exactly 2. In the generic case of irrational values of a, Kesten's 
question leads to an interesting number-theoretic problem. Some results can be found in his 
original paper |6j where it is stated that for almost every a positive measure set of orbits escapes 
to infinity but return to zero infinitely often. Most contemporary analysis of this phenomena 
can be found in |llj . 

Surprisingly, Kesten's system becomes important in the study of discontinuous twist maps 
after an appropriate renormalization procedure is carried out. 

Elementary properties of az-map. For non-degenerate twist z 7^ the following properties hold: 

• For any z < — 1 nearly half of trajectories of the system ^ escapes to infinity. It 
immediately follows from the fact that converges. 

n 

• Fix z E N then for a = 1, it is easy to verify that half the orbits are unbounded and for 
a = ^ all the trajectories are periodic. 

• The most interesting and difficult problem of boundedness occurs for z G [—1, 1). 
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3. PiNBALL SYSTEM 

Consider now Fermi- Ulam like system with the fixed walls, but with one of the walls containing 
a pinball mechanism: the momentum of the particle increases or decreases when it hits the wall 
according to the following law: 



(3) 




v + liite [0,1/2] (mod 1) 
v-liit£ [1/2,1] (mod 1), 



i.e. the momentum is increased (decreased) during the first (second) half period. This dynamics 
is described by the map with a being the distance between the walls and z = — 1. 

We rewrite the system ^ for z = — 1 and we will call such system a pinball transformation 
that will be denoted by P. For the sake of clarity, it would be more convenient to consider the 
base circle (p G [0, 2). 

(4) /^i = (^+t) "'"^^ 

Ui =/ + sgn((/Ji -1) 

Numerical experiments show that for typical values of parameter, the trajectory of the sys- 
tem is nearly recurrent for a long time and moreover approximates some piecewise smooth 
function having singularities only at the discontinuity lines ip = \ and 99 = 0. 




Figure 1. Left: typical trajectory in a phase space: 10 iterations starting at 
I = 50. Right: Fundamental domain <1>. 

Our goal is to explain this behavior by renormalizing the induced transformation in the so- 
called fundamental domain The fundamental domain is related to Poincare' section for flows 
in the sense that any orbit returns to it. In the Pinball map, we define the fundamental domain 
as the set of points ((^, /) laying between singular line 99 = and its image 

(5) ^. = {((^, /) G [0, 2) X [0, 00], < < Proj^P(0, /)}. 

Our first statement concerns the asymptotic description of the first return map: 

Theorem 1. First return map of the domain $ under the map ^ is a 0{I^^) perturbation of 
the Kesten transformation ^ with rotation number // = expa~^, provided a > l/ln3. 

If ^ is rational, than the reduced Kesten system might possesses a uniformly growing trajec- 
tory. Indeed, if one of the periodic points stays longer in the positive part of the base interval 
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than in the negative part, then the corresponding trajectory grows without bound. Our con- 
struction of an escaping trajectory of the system Q consists in choosing the initial data in an 
appropriate way so as to kill the leading order perturbation of Kesten's map. Next, we would 
have to estimate that the remaining perturbation will not destroy such "resonant" growth. 
Combining these ideas, we prove 

Theorem 2. For a = — ^ there exists an unbounded trajectory in the system Q. 

Note that smoothing the signum function discontinuity in Q will make KAM theory appli- 
cable and then all solutions will be bounded. 

Indeed, change the variables: {^,y) = {ip,I~^). In the new variables, the smooth version of 
the Pinball transformation takes the form 

= ip + ay mod 2 
1 



where / is smooth and \ f{<p,y)\ < 1. Then 

1 

^'^-yf{f,y) ^ 



l + y(yf{^,y)y 



n=l 



SO the perturbation is of order 0{y'^) which is much smaller than the twist. The curve intersection 
property follows from the area-conservation in the original variables. Therefore, this map satisfies 
the conditions of monotone twist theorem, see e.g. p.Oj. 

4. Proof of Theorem [H 
Define a subset <I> C S x M_|_ between the discontinuity line (f) = and its first iteration: 

(6) <i> = |(^,I)|(/.G fo, " 



- 1 

and consider the transformation T({p, I) = {ip' , I') as the first return map for any point {if, I) G ^ 
according to Q. We have the following bound on the action change 

Lemma 1. If(ip',I') is an image of the point {(p, I) under the transformation T then \I' — 1\ ^ 1. 

In other words. Lemma [T] states that as the angle variable winds around the cylinder and the 
action variable undergoes large changes, after returning to the fundamental domain, the action 
will not change by more than 1. This statement gives a good local control on the orbits. 




Figure 2. Mapping from fundamental domain to itself represented on the base 
circle (right panel) and on the covering space (right panel). 
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Next we describe the structure of the subsets in the base (p E (0, a/(/ — 1)) for which the 
action increases I' = I + 1 or decreases I' = I — 1. 

Let us introduce the resettled angle variable (p = {I — l)(/9a~^. In the renormalized variables, 
the fundamental domain can be represented by 

$ = ((^,/) = [0,l]xM+. 

Lemma 2. The set <I> n {/ = C} is the union of three subsets 

$n{/ = C} =X+UX_ UXo, 

where and I- ttre intervttls of equal mettsure and consist of all points {(p, C) such that 
T((/9, C) = (•, C lb 1) respectively. 

The intervals are contained in the regions 

^ I a \ , ^ f a(l — n~^) a 

X+ e 0, — — - and X_ e ' ^ ^ ' 



fx{I-l)J V ^-1 '^-1. 

where ii{a) = exp(a^"'^). 

Finally, we have the following lemma which ends the proof of Theorem [T} 
Lemma 3. //a > l/ln3, then the first return map T in the resettled vttrittbles tttkes the form 

>i = ^ + 2/^(a) + 0(/-i) 
h = I + X+i'Pi) - X-i'fi) 



Tip, I) 



where //(a) = exp(a x± o.^^ chttrttcteristic functions of positive and negative intervals de- 
pending weakly on I. 



5. Proofs 

5.1. Proof of Lemma [l]. 

Proof. We begin with giving a heuristic argument based on the Figure [2] (right panel) . We can 
represent orbits as stairs going up and down in the domain (p S [0,2]. Consider the special 
configuration that is symmetric with respect to p = 1. It is easy to see that the corresponding 
solution will have the same action after the first return map. By moving the graph so that the 
top level does not cross = 1 it is possible to see that the change in the action cannot be more 
than one. Indeed, the lowest level steps are wider than the top one and therefore at most one 
crossing can occur. 

Now, we provide the full proof. Assume the initial point (p, I) is in the fundamental domain 
i.e. p S (0,a/(/ — 1)). The map Q is iterated n + 1 times while / increases until the orbit 
is one step away from crossing p = 1. Next, the map is iterated n' + 1 times while / decreases 
until the orbit is one step away from crossing p = 2. Then, by the definition of the fundamental 
domain, last (n + n' + 2)-nd iterate is in $ and we obtain 

(7) T{p, I)=(p + aS,{p, I) + aS2{p, I) + + " , .. J + n - n') , 

\ l+n+ll + [n — n' — 1) ) 

where 

n ^ n' ^ 

k=0 k=0 
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The numbers n and n' are uniquely defined by the relations: 

fl > ip + aSi{ip,I) 



(9) 
and 

(10) 



|l < <f + aSi{ip,I) + 

2> ip + aSi{ipJ) + aS2{ip,I) + 
2< ip + aSi{ip,I) + aS2i>f,I) + 



a 



I + n + l 



a 



I + n + l 

a 



+ 



a 



I + n + 1 I + n — n' — 1 

The equation Q means that n-th iterate is the last one staying in the right half-circle, so the 
next iterate will be in the left half-circle. Similarly (n -|- 1) -|- n'-th iterate is the last one before 
returning to see Figure [2] 

Denote by AS" = 52 — 5*1. Rewriting the second sum in ^ for k' = n-k we obtain for A5: 



(11) 



A5 



n ^ n ^ 

^ I + k' I + k' 



k'=n~n' 



Note that the expression (11) implies that 



(12) 



A O 1 1 

A5 > r + 



A5 < 



I-l 1-2 
1 1 



fc=0 

if n — n' < —1 
if n — n' > 1. 



/ / + 1 

To prove that |n — n'| ^ 1 it is sufficient to show that (12) contradicts ^ - (10). Rewriting 
([lOj) for 52 = A5 Si we obtain 

CK O. 

(13) 2 >(^ + 2a5i(v9,/) + aA5+ ^ ->2 



I + n + l 



I + n — n' — \ 



Multiplying the equation ([9|) by 2 and substituting into (13) we get 

1 



(14) 



AS + 
AS > 



a I + n + 1 

V9 1 



1 



a 7 + ^^+1 I + n — n' — 1 



Since for any point (ip, I) £ ^ one has — G ( 0, 

a V I — 1 



5.2. Proof of lemma [21 



inequalities (14) contradict (12). 



□ 



Proof. As already found in the proof of Lemma [T| transformation T in 99- variable takes 

the form 



(15) 
where 

(16) 



^' = ^+ 2a5i(^,/) + 



a 



I + n + l 



2 + aA5, 



A5 



1 



if 99 G 
if 99 G Xo 



I-V 
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Define 6i, 62 by the relations (see Fig. [2]) 

(f + aSi{(p,I) = 1 - 61 

o. 

+ aSi{ip, I) + — — = 1 + 62, 

i + n + 1 

then it is easy to see that 

Lp + 2aSi{ip,I) + 



a 



/ + 1 



5i - if. 



Then the expression ( 16 ) could be rewritten in the form 

for G X+ 

-,0^ , for (/? G Zq 

for if G X_ 





'62 


-Si 


-ip>0, 


(17) 


52 


-Si 


- G (- 




S2 


-Si 


if < 



a 

a 



I -I' 



In particular, the first inequality in (17) implies that for 93 G Xj^ one has 

a 



^ < 



a 



smce 



I + n + 1 

a 



S2 + S1. 

a 



I + n + 1 

Similarly, it is easy to verify that if G X_ then 09 > 

I — 1 / + n + 1 

describes some rigidity properties of the intervals Iq , X_|- , X_ in (|17[) 



. The next proposition 



Proposition 1. Assume that {ipQ,I) G X+ and that the corresponding sum Si{ip(),I) has n 
terms. Then the point [ip^ = ipo + e, I) belongs to X+ if and only if the sum Si{ipe, I) has also n 
terms. 

Proof. Denote by S^"^ and S^"^ the corresponding parts of the interval ;— ;- for the point 



(/7o- By assumption lpq G X_(_ and thus S^'' — S^"* — (/?o > 0, then by (17) it is sufficient to verify 

that S2^ — S^'^ — i^pe remains positive. 

Then the proof can be obtained from the following calculation: 

(a) First, consider the case when the number of terms remains the same (equal to n). Then 

S2^ = 82 + e, S^^ = 5i — e and so 



I + n + l 



(18) - jj^) - ^ = (4°) + e)- (4°) -e)-e = - sf^ + e > - sf'^ 

(b) Assume now that the number of terms in Si{ipe, I) is different from n. Assume it contains 
1 terms (all other cases can be treated similarly). Then we have, see Figure [s} 

1 



n 



I+n+1 




Figure 3. Rotation by e. 
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I + n + l 
which imphes 

= ^-4^) = ^-6+4°). 

^ I+n ^ I+n ^ 



Finally, we have 



I+n+l I+n 
The latter expression is negative since by construction 



^ I+n I+n+l 

This ends the proof of the proposition. □ 

Using proposition [T| we will find the set X+ . Let the initial angle be 999 = then 



1 1 




n+l 

4°) = laY-^) and 4°) = -^^ sf^ 



If 4*^^ > 4*^^ then by (17) the point (0,/) belongs to and then proposition [l] implies that 

x+ = (o,4°0- 

Otherwise, if 4'^^ < ^1'^^ then (0, /) does not belong to X+ and to find the leftmost point 
{4> = G 2^+ we need to satisfy the condition 4*^^ ~ '^1'^^ — e > 0. Using (18) we obtain 
e = — 4*^^ and therefore in this case X+ = isf'^ — 4*^^ 1 4*^^ 



Similar calculations can be carried out for X_ , but it is easier to use symmetric properties of 
the map with respect to reversing the "time" . 

Then, if (5f °^ > 4°°^ (4°°^ correspond to but are obtained from negative iterates) the 
point (0—,/ — 1) which is the preimage of {j^ — ,I) belongs to X_ and by the same argument 
as above for X4., we have 



I-l ""^ '/-I 



If the converse holds, i.e. 6^^^ > 6f'^^ we obtain 



This ends the proof of the lemma. □ 



The location and the measure of the intervals and X_ are then controlled by the fractional 
parts of the sums a5i(0, /). We shall estimate these quantities in the proof of the next lemma. 

5.3. Proof of Lemma [si 

Proof. For the clarity of presentation, we give the proof only for the case I = N. The general 
case can be treated similarly. We introduce the notation for the expression 



aV(^,/) = 2a5i(v9,/) + 

I+n + l 



from (15). 
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In the rescaled variables, the transformation (15) takes the form 

1 



(19) 



a r 
a r 
rv I 



1 + 



If + I'ip{I,ip), 



I -I 



1 



1 



1 



<p+{I-2)^P{I,ip)+2 



I -I 



(pel- 



Case if £ X+. If the action variable / increases, then the number of iterates of the transformation 
Q contained in the right half-circle (where the action is gained) is greater than the number 
of iterates contained in the left half-circle. Lemma [T] implies that n = n' + 1. Thus, the 
transformation T for these values of ip takes the form 



ip + 2aSi{ip,I) + 



a 



I + n + l 



where 



5i((^,iV) = ^ 



1 



fc=0 



N + k 



H 



N+n 



H 



N-l, 



where denotes k-th harmonic number. The harmonic number has the following asymptotic 
expansion 



Hk=lnk + j + 



2k 



2j 



where Bj denotes j-th Bernoulli number. Thus, 



(20) 



, 7V + n 1 

: In h - 

iV- 1 2 



1 



2jfc2i 



1 



N - 1 
+ 



1 

iV4 



N + n 

12 \{N + n)^ {N-iy, 

We will now show that n = [fJ,{N — 1)] — A^, where fi = exp(l/a). By proposition [T| the 
number of steps in the positive half of the cylinder n is independent of (/? G X_|_, therefore it is 
sufficient to verify n = [iJi{N — 1)] — N only for some (p using ([9]). Moreover, we can take 99 = 
even though this point might not be in X+. Indeed, according to lemma [2] the point ip = 
satisfies p € X+ or G Xq. If the former, we are done. If the latter then it is easy to see that 
by increasing (p we will eventually cross into X-|- without changing n: the number of steps in 
positive part of the cylinder. 

We have 



iV- 1 



1 



+ 



1 



2{N 
( (N- 



1) 
I? 



N - 



1 



1 + 



1 +0 



1 

w 



Using the relation fi{N — 1) 



12(iV-l)2 V[/u(Af-l)]2 
[n{N - 1)] + {fi{N - 1)} and denoting y : 
1 / 1 



rewrite 



1 

12(iV- 1)2 V/i2(l-y)2 
Expanding in perturbation series for y we obtain 



Mi-y) 

1^0 



- 1 
1 

W 
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1 



(21) 



1 



k=l 



2{N -1) \n 
1 



12(iV- 1)2 



1 + 



2/x(iV-l) 12;u2(iV-l)2 k 



+ 



Recall that a = and 1 < /i < 3. Then using (21) it easy to verify that aSi{(p,N) < 1 



in the limit N ^ co as the leading order terms containing are negative. On the other 
hand, J^_^}^_^l = [^(Af-i)]+i ^^"^ similar estimates aSi{ip,N) + > 1. This shows that 



Af+n+l [ti{N-l)] 

n = [^i{N — 1)] — is the correct number of steps. 



Now, we will derive an explicit expression for the first return map T. Multiply (21) by 2a 

1 



(22) 



2aSi{^,N) = 2 + 



k=l 



a 



1 



a 



1 

6(iV- 1)2 V/? 



(N-l) V/x 
1 

fiiN-1) ~ 6/x2(Ar- 1)2 ~ k 



1 + 



1 



1 



+ 



1 

W 



so the transformation T takes the form 



92 



ip + 2aSi{ip,N) + 



a 



(23) 



9^ + 



a 



(N-l) 



1 



1 



- 1)] + 1 



a 



mod 2 



fe=i 



1 



1 

6(iV- 1)2 V/? 
1 



1 + 



/i(Af-l) 6^2(^_i)2 A;)^'^(iV4 



l + [MiV-l)] 
mod 2. 



+ 



. — . N-l 
Multiplying both sides of the equation (23) by the scaling factor and using the expansion 



a 



N-l 



(l + [M(iV-l)]) /xll 



we obtain 



N-l , 

-f = f + 



a 



(24) 



00 / 



1 



y 



1 



1 °° / 



k=0 



1 



lj{N - 1) + 1 



fc+i 



k=l 



6(iV- 1) V/W^ 



fi{N -1) + 1 



+ 



Mi + Miv-1)) 
1 1 



+ 



+ 



H{N-1) 6^2(jY-l)2 k 



Since (/? G X+ we can use ( 19 ) 



1 \ iV- 1 , 
1 + I V'- 



a \ N — 1 J a 
Thus, in the leading order in the rescaled variables, the transformation takes the form 

'2 



1 . 



Remark 1. Note that in case (/? G Xq or ip £ it is not necessarily true thatn = [n{N—l)]—N, 
however the error of such approximation e{N,fj,) = |n = [^{N — 1)] — N\ is bounded by a constant 
independent of N. Such small correction does not have an effect on the leading order asymptotics 
below. 
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Case G Xq. In this case n = n' and so rescaling for the image (f' coincides with the rescahng 

for ip. On the other hand in this case AS = so, for the rescaled variables we get from 

2 



Q : (^' = (^ + - +0(iV-i). 

Case 09 G T_. In the third region n = n' — 1 and so AS* = \ . Rescaling 

gives id' = if' = I 1 I tp' + 2 H which again gives from (24) (5' = 

a V N — 1 J a N — 2 ^ — ' 



(^ + - + iV-i . 

□ 

5.4. Proof of Theorem m 

Proof. If a = then jj, = n and therefore y = 0. Thus, all the terms in series expansion (24) 

Inn ' — ' 
vanish. In the leading order, the rescaled transformation takes the form 

= (- -1 

i.e. the rigid rotation on the circle. 

For our particular choice /i = 2, in the leading order the transformation T is the identity 
mapj^ Thus, we get 

To construct an unbounded orbit, we search for an initial point ifo satisfying two conditions: 

• X+i'fo) = 1, 

• The image <^q under the transformation T coincides with initial angle up to the new 
normalization (fQ. 

Note that for a larger action variable the scaling factor would be different. Therefore, we have 

N-1 1 

and the "renormalized fixed point" condition yields 

- = + 8(iV^ - 4iV^ + 
Thus, up to the order 0{N^'^) 

(fo AN -2-8N + 8 2N -3 



N 8{N-l){AN-2) 8{N-1){2N-1) 

1 11 1 



8(A^-1) 4(7V- l)(2iV- 1) 8A^ 8N{N-1) A{N - 1){2N - 1) 
1 1/1 1 \ 1 1 



8A^ 4(iV- 1) V2iV 2iV-iy 8N 8N{N - 1){2N - 1)' 
As a result, for (po = - + 0{N~'^) its image under first return map is given by (p^ = — ^0{N' 



^For other values /i £ N, we would have periodic mappings. 
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To check the first statement x+iv^o) = 1 it is sufficient to show that ipo + 2aSi 



a 



is greater than 2. This immediately fohows from (22) and the estimate 

1 



N + n+l 



ipo + 2aSi + 



a 



> 



N + n + 1 8iV hi 2 



+ 2 



1 



hi2(27V-2; 



+ 0{N-') + 



1 



{2N - l)ln2 



> 2. 



□ 



iJllplK" 



■ ■ 



Figure 4. The structure of the fundamental domain for a 



1 



log 2 



. The regions 



of positive growth are in red and their images are in green. The regions of 
negative growth are in blue. 



Remark 2. Construction of the unbounded orbit in the Theorem^can be easily extended for any 
a = i^^2m ' £ I^- iemmasJT] an(i[£| do not depend on a and can be applied without any modifi- 
cations. The only difference with the considered case ^ = 2 is in the action of T on the neutral 
setlf). Instead of the pure rotation, there will be an additional interval exchange transformation 
acting on Iq with [fj] — 2 subintervals. We omit technical details of this construction. 

Remark 3. We conjecture that the same arguments as in Lemmas\^ -[^ can be used for other 
values of z € (—1,0). Here instead of asymptotic expansion for harmonic numbers one should 
use generalized harmonic numbers Hn,z- Since Hn,z C(^) as n ^ oo one can establish the 
relation for -ff(/+n),z ~ -f^(/-i),z = 2 to find the expression for n. Arguments in Lemma^ and 
can be applied with minor changes. 

Appendix A. 

In the appendix we give more detailed derivation of some specific problems that lead to az 
map. 

Example 1: Particle in switching potential. Consider a classical particle moving on the 
line in the square wave potential V{x) = (— and assume the potential is switched every 
time unit V{x,t) = (—1)'*^ • (— l)'^'. While such potential is not differentiable, there is a natural 
way to define the dynamics by using the energy relation: the kinetic energy changes by 2 if 
the particle passes t £ Z integer points. We should ignore the singular subset of the extended 
phase space {t, x, x) where there is discontinuity in both time and space and the dynamics is not 
defined. Such subset has zero measure. Outside the singular set particle moves with constant 
speed v{x,t) = \JE — V(x^i). 
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Since the dynamics defined for x mod 2, t mod 2 we can project the dynamics on a plane 
onto the system on a cyhnder [0,2) x M+. Write the Hamiltonian form of the unit time step 
transformation for this system. Hamiltonian has a form H{x^ x, = ^ + V'{x^ t) and so for the 
canonical action-angle variables we get 

•l± = ^2{H - I) + y'2{H + 1) 
(25) <i ^ ^ / x-1 1 \ 

'^^ T l^v^2(i7±sgn(rE-l)) V2(^ T 1) 7 

H — I is a period of rotation along the level curve H(x, x) = H 



Here T = — , , , 

\^2{H-l) ^2{H + l) 

and {'^±-,I±) correspond to the angle-action variables for odd/even values of [t]. To deduce the 

system in action-angle variables, write down a generating function. For Hamiltonian we have 



Where 



Or, using (25) 



H{x,x,t) = H^{I,t) + d^S±{x{t),I,t) 



S+{x,I) = {x - 1)^2(1 + sgn(x - 1)) + ^2{I - 1) 
S-{I) = iS+{-I) 



2r(/+l)<^+ 



72(7^ 



:, X > 1 



_2r(/- i)v9+, X < 1 



Finally, the total system in variables {(p,I) can be deduced from the expression (25) and the 
relation ip± = djS±. 



(26) 



(fi = {(f + aVT) mod 2 
h = / + sgn((^i - 1) 



where a is some constant. Clearly, this system can be considered as a particular example of 
transformation ^ with z = ^. 

Example 2: Fermi-Ulam acceleration. In Fermi-Ulam problem the particle bounces be- 
tween two walls. Assume that one wall is at rest x = and the other moves periodically 
X = p{t), p{t + 1) = p{t) > 0. There is a standard transformation "stopping" the wall, see e.g. 
[Ti] with 

ds ^ 

x=p{t)y, T = I -Yf-rds. 

Jo P [s) 

In the new variables, the equation takes the form 

" I •■ 3 n 

y +PP y = 0, 

where / denotes the derivative wrt r. Evaluating the one period map, under the assumption 
that p is piecewise linear, we obtain 



(27) 



y2 = {yi + y'i) modi 

y'2 = y'l +y2Sgn(y2 - 1/2). 



This mapping is not a particular case of az-map but it corresponds to the linear growth of 
the action when z = 1. While showing unbounded growth is relatively easy in this case, more 
challenging problem is to estimate the relative measure of bounded solutions. This has been 
done in f2j. 
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Example 3: Outer billiards with degenerate boundary. Consider the outer billiard sys- 
tem. Let 7 be a smooth strictly convex curve on the plane. Take any point x G outside of 
7 and let l{x) be a ray tangent to 7 and oriented in the counter-clockwise direction. There is 
another point on the ray T{x) which has the same distance to the tangency point as x. This 
defines the outer billiard map, see e.g. |13j . A natural question is whether all the orbits are 
bounded. Thus, one is led to study this map for large x. 

Assume that 7 is a unit circle centered at the origin, then for large x the square of the map 
is close to identity and it leaves concentric circles invariant. The angle changes by a factor of 
l/l^l, which corresponds to z = — 1 in the az-map. Now, consider the circle with a small circle 
segment removed. Then, the map becomes a small discontinuous perturbation of the above 
integrable map. When, 7 is half the circle, the outer billiard has unbounded orbit, see [HIS]. 

Earlier, Schwartz [T2j constructed unbounded orbits for quadrilateral 7 . 
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